CAN MERONS DESCRIBE CONFINEMENT? 



JAMES V. STEELE 
Center for Theoretical Physics and Laboratory for Nuclear Science 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA 
E-mail: jsteele@mit.edu 

Merons, conjectured as a semiclassical mechanism for color confinement in QCD, 
are topological charge-1/2, singular solutions to the classical Yang-Mills equations 
of motion. I will discuss how lattice techniques can extend the study of merons 
to nonsingular stationary solutions without destroying properties believed to be 
essential for confinement, and how zero modes can be used to identify these gauge 
field configurations in stochastic evaluations of the lattice QCD path integral. 

1 Introduction 

Merons am one of the earliest proposed semiclassical mechanisms foL con- 
finement, El but have received less attention than other approaches Ha due 
to analytical limitations. I will present a new direction for meron studies ,0 
developed recently, □ in which lattice gauge theory will be able to determine 
the presence and role of merons in QCD. 

Chiral symmetry breaking and confinement of color charge, the two essen- 
tial features of low-energy QCD, cannot be understood perturbatively in the 
coupling constant. Two complementary nonperturbative approaches are an- 
alytical approximation using semiclassical analysis of the QCD path integral 
(analogous to the WKB approximation in quantum mechanics) and numerical 
solution using lattice gauge theory. Combining physical understanding from 
the semiclassical approximation with systematic rigor of the lattice has pro- 
vided valuable insight into hadron structure B and can be extended to study 
confinement. 

Most semiclassical treatments of QCD have focused on instantons, leading 
to a qualitative, and in some cases even quantitative, understanding of chiral 
symmetry breaking and the structure of the QCD vacuum. □ However, simple 
configurations of instantons, such as a dilute gas EJ or a random superposi- 
tion, Q do not confine color charges. It is possible the complicated dynamics 
of confinement escapes any explanation save a full solution of QCD at large 
coupling. However, analogy with lower dimensional models seems to indicate 
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otherwise. 

In 2+1 dimensions, the Georgi-Glashow model, which couples a triplet 
Higgs to an SU(2) vector boson field, contains magnetic monopoles as solutions 
to the Euclidean classical equations of motion. Taking the Higgs to have an 
expectation value along the third isospin direction (jia — <5a3 v, the residual 
massless photon exhibits a monopole solution of form ^ £3^1^ x^jx^ at large 
distances. Polyakov showed these monopoles condense, leading to confinement 
of electric charge. 13 Similarly, the Schwinger model (QED in 1+1 dimensions) 
exhibits confinement through the dominance of vortex configurations □ of 
form — t^iyXvjx^ . These vortices have a topological chaLge of i, which 
also shows up in the quark correlator at finite temperature. l2I 

These semiclassical examples of confinement are intriguing in that the ef- 
fect is generated by the same gauge-field configuration interpreted in different 
dimensions. The meron is the (3+l)-dimensional nonabelian extension of the 
't Hooft-Polyakov monopole configuration with = rja^iuXu / x^ . The impor- 
tant property shared by each of these configurations is A ^ 1/r and is not 
pure gauge. Then the integral §A^dx^^ over a loop of area R x T will have a 
contribution of order unity as long as the source of the field is within R of this 
loop in any direction. El If the vacuum is made up of a uniform distribution of 
these semiclassical configurations, the energy between charged particles will 
behave like E{R) oc i?'* in d+1 dimensions. Interactions between monopoles 
in 2+1 dimensions reduce this down to a linearly rising potential u and it is 
conjectured that merons will do the same in 3+1 dimensions. EJ 

However, a purely analytic study of merons has proven intractable for 
several reasons. Unlike instantons, no exact solution exists for more than two 
merons, llil gauge fields for isolated merons fall off too slowly (A ^ l/x) to 
superpose them, and the field strength is singular. A patched Ansatz configu- 
ration that removes the singularities 111 does not satisfy the classical Yang-Mills 
equations, preventing even calculation of gaussian fluctuations around a meron 
pair. E3 

The properties of merons are reviewed in Sec. 2. The discussion is re- 
stricted to SU(2) color, since topological objects live in embedded SU(2) sub- 
groups of larger gauge groups. As work with instantons has shown, fermionic 
zero modes are a powerful tool for identifying topological objects on the lat- 
tice. Therefore, in Sec. 3, the zero mode of the continuum meron pair Ansatz 
is analytically derived. In Sec. 4, a smooth, finite action configuration corre- 
sponding to a meron pair is presented, which with its associated quantum fluc- 
tuations, should be present in Monte Carlo lattice QCD calculations. These 
exact numerical results along with the associated zero mode are compared 
with analytical results of the meron pair Ansatz. 
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2 Analytic Meron Review 



Two known solutions to the classical Yang-Mills equations in four Euclidean 
dimensions are instantons and merons. Both have topological charge, can be 
interpreted as tunneling solutions, and can be written in the general form (for 
the covariant derivative = d^ — iA'ta'^ /2) 



(1) 



with f{x^) = lipP' + p^) for an instanton and /(x^) = i for a meron. 

Conformal symmetry of the classical Yang-Mills action, in particular un- 
der inversion — > x^/x'^, shows that in addition to a meron at the origin, 
there is a second meron at infinity. These two merons can be mapped to 
arbitrary positions, defined to be the origin and d^, by the conformal trans- 
formation!^ 



,2 - C?)^ 

(X - dY 



(2) 



After a gauge, transformation, the gauge field for the two merons takes the 
simple form 



x"" {x - df 

X^ 



(x - dy 



(3) 



Similar to instantons, lI a meron pair can be expressed in singular gauge 
performing a large gauge transformation about the mid-point of the pair. 



resulting in a gauge field that falls off faster at large distances (A 



A 



camful treatment of the singularities shows that the topological charge density 



Q{x) 



Tr 
16^ 



F,.F^'^]=\s\x) + \s\x-d) 



(4) 



yielding total topological charge Q ^ 1, just like the instanton. 

The gauge field Eq. (0) has infinite action density at the singularities 



''A more general conformal transformation only differs from this choice by either a rotation, 
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= {0, dp}. Hence, a finite action Ansatz has been suggested 







a;^ < r , 



region I, 



1 

2 



vo/ii/-^ , r < Vx'^ < R , region II, 

R < Va? , region III, 



(5) 



I a;2 + i?2 ' 



with arbitrary radii r and R. Here, the singular moron fields for \l x^ <r and 
\fx^ > R are replaced by instanton caps, each containing topological charge 



to agree with Eq. (^. The action for this configuration is 

o = — TT- TT- In — , 

g 9 r 



(6) 



which shows the divergence in the r — > or i? — s- cx) limit. There is no 
angular dependence in this patching, and so the conformal symmetry of the 
meron pair is retained. For example, under a dilatation a;^ Xx^, both 
r and R get multiplied by 1/A but the ratio and hence the action Eq. 
remains invariant. Although this patching of instanton caps is continuous, 
the derivatives are not, and so the equations of motion are violated at the 
boundaries of the regions in Eq. (H) . 

Applying the same transformations used to attain Eq. (^) to the instanton 
cap solution, regions Land III become four-dimensional spheres each contain- 
ing half an instanton. □ The geometry of the instanton caps is shown in Fig. |l] 
for the symmetric choice of displacement d = \/ Rr along the z-direction. Note 
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that the action for the comphcated geometry of Fig. |] is still given by Eq. , 
which for d ~ \J Rr can be rewritten as 

^ = 5o(l + Jln^) , (7) 

with — Sir'^/g'^. This case will be used below, and generalization to a 
different d is straightforward. 

The original positions of the two merons = {0, d^} are not the centers 
of the spheres, nor are they the positions of maximum action density, which 
occurs within the spheres with ^max — (48/g^)(i? + r^/d^ at 

(^^i). = , (A/„i)^ ^ d, . (8) 

However, in the limit r ^ and i? — s- oo holding d fixed, the spheres will 
shrink around the original points reducing to the bare meron pair in Eq. (^. 
In the opposite limit R — > r, the radii of the spheres increase to infinity, 
leaving an instanton of size p — d. 

An instanton can therefore be interpreted as consisting of a meron pair. 
This complements the fact that instantons and antiinstantons have dipole 
interactions between each other. If there exists a regime in QCD where meron 
entropy contributes more to the free energy than the logarithmic potential 
between pairs, like the Kosterlitz-Thouless phase transition, instantons will 
break apart into meron pairs, u with the intrinsic size of the instanton caps 
determined by the original instanton scale p. This will occur when the entropy 
to create meron pairs, which is R/r in each space-time direction with the 
general scales r and R given in Eq. (|^) , overcomes the logarithmic interaction 
energy Eq. (||) to dominate the QCD path integral, or schematically. El 

How the scale s in the renormalized coupling constant depends on r and R is 
not known. It is assumed that as longer and longer distances are probed, s also 
increases, causing entropy to win out over energy in Eq. (^. This occurs for 
relatively weak coupling, f/^(s)/87r^ > so the semiclassical approximation 
should still be valid. El Lattice calculations can test this physical picture by 
increasing the coupling constant and seeing if the predominant semiclassical 
degrees of freedom change from a dilute gas of instantons into a plasma of 
merons. 
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3 Analytic Zero Modes 



The Atiyah-Singer index theorem states that a fermion in the presence of a 
gauge field with topological charge Q, described by the equation 

pip = Xip, with ip=(^^^ , (10) 

has rifl right-handed and ul left-handed zero modes (defined by A = 0) such 
that CL — TiL — nji. This can sometimes be strengthened by a vanishing theo- 
rem. E3 Applying twice to ip decouples the right- and left-hand components, 
and focusing on the equation for ipn, gives 

+ ^Vaf^.T'^F^''^ ^R = 0, (11) 

where fja^i, denotes the 't Hooft symbol, cr" acts on the spin indices of ipRi and 
F^'^ acts on the color indices. For a self-dual gauge field (like an instanton), 
the second term in Eq. ([ll] ) is zero; and since is a negative definite oper- 
ator, there are no normalizable right-handed zero modes, implying Q — n^. 
Although the meron pair is not self-dual, the second term can be shown to be 
negative definite as well, leading to the same conclusion. 

In general, a gauge field in Lorentz gauge with Q = I can be written in 
the form 

^m(^) = -^w-^'^i^n(x) . (12) 

This has a fermion zero mode given by 

, C-AAn-V^e^ (13) 

with normalization J\f and e = ia2 coupling the spin index a to the color 
index a (both of which can be either 1 or 2) in a singlet configuration. The 
gauge field for the meron pair with instanton caps can be expressed in the 
form Eq. ( |T2| ) with 

n(-) - <! ; (14) 

region II, 



■\/a;2(x — d)^ 
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with - 
then Eq. 



r/d and 
3h with 



R/d. The normahzed solution to the zero mode is 



2 - , — 
R 



1/2 



(15) 



Note that the unxegulated meron pair (r 
zero mode itself E3 



del 



2\/27r(x2(a;-d)2)3/4 



0, i? oo) has a normalizable 



(16) 



The gauge-invariant zero mode density ilj''ip(x) has a bridge between two 



merons that falls off like 



-3 ; 



in contrast to the 



fall-off in all other di- 



rections. This behavior can be used to identify merons when analyzing their 
zero modes on the lattice, similar to what was done for instantons in Ref. [16]. 



4 Merons on the Lattice 

As mentioned above, the patching of instanton caps to obtain an explicit 
analytic solution has unavoidable and unphysical discontinuities in the action 
density. Therefore, in order to study this solution further, the gauge field is 
put on a space-time lattice of spacing a in a box of size Lg x ii x L2 x L3. The 
gauge-field degrees of freedom are replaced by the usual parallel transport link 
variable. 



C/^(x) = Pexp 



A^{z)dz'' 



(17) 



The exponentiated integral can be performed analytically within the instanton 
caps, producing arctangents. Outside of the caps, the integral is evaluated 
numerically by dividing each link into as many sub- links as necessary to reduce 
the O(a^) path-ordering errors below machine precision. 

Calculating the action density S{x) for a meron pair with instanton caps 
using the improved action of Ref. [17] and the topological charge density Q{x) 
using products of clovers, the lattice results (open circles) are compared with 
the patched Ansatz results (solid line) in Fig. |^ for r = 9 and d = 20.8 (in 
units of the lattice spacing) on a 32'^ x 40 lattice. The Arnoldi method is used 
to solve for the zero mode of this gauge configuration and hence the density 
'0^i/'(^)i which is also compared with the analytic result in the same figure, 
showing excellent agreement. 

Two important features of the patched Ansatz also evident with the lattice 
representation are the discontinuities in the action density at the boundary 
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2e-04 




z 



Figure 2. The action density S{x) for the regulated meron pair, normaUzed by the instan- 
ton action So = Sn'^/g'^, sliced through the center of the configuration in the direction of 
separation. Shown are the analytic (solid), initial lattice (open circles) and cooled lattice 
(filled circles) configurations. The same is shown for the topological charge density Q{x) 
and fermion zero mode density ip'^^{x). 



of the instanton caps and vanishing of the topological charge density outside 
of the caps as given by Eq. (Q). The discontinuity in the action density is 
unphysical and can be smoothed opi by using a relaxation algorithm to itera- 
tively minimize the lattice action. E3 On a sweep through the lattice, referred 
to as a cooling step, each link is chosen to locally minimize the action density. 
Since a single instanton is already a minimum of the action, this algorithm 
would leave the instanton unchanged (for a suitably improved lattice action). 
The result for the regulated meron pair after ten cooling steps is represented 
in Fig. 1^ by the filled circles, showing the discontinuities in the action density 
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Figure 3. The topological charge density Q{x) on the left and formion zero mode density 
tp'^tp{x) on the right for a cooled meron pair with r = 9, separated by distance d = 20.8, 
17.7, and 14.1 respectively, in the (z,t)-plano. 



have already been smoothed out. Figure^ shows side- by-side the cooled Q{x) 
and ip^ipix) in the {z, f)-plane for three different meron pair separations. As 
the separation vanishes, the zero mode goes over into the well-known instanton 
zero mode result. 

Like instanton-antiinstanton pairs, however, a meron pair is not a strict 
minimum of the action, since it has a weak attractive interaction Eq. and 
under repeated relaxation will coalesce to form an instanton. This is analo^ 
gous to the annihilation of an instanton-antiinstanton pair through cooling. El 
Nevertheless, just as it is important to sum all the quasi-stationary instanton- 
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antiinstanton configurations to obtain essential nonperturbative physics, 
meron pairs may be expected to play an analogous role. A precise framework 
for including quasi-stationary meron pairs is to introduce a constraint Q[A] 
on a suitable collective variable q (here chosen to be the quadrupolc moment 
3z^ ~ ~ ~ of the topological charge) as follows 

Z = JdA cxp{-S[A]} = jdqZq, (18) 

= jo A exp [-S[A] - A [Q[A] - qf^ . (19) 

The meron pair is then a true minimum of the effective action with constraint, 
allowing for a semiclassical treatment of Zq. Afterwards, q is integrated to 
obtain the full partition function Z . 

The criterion for a good collective variable q of the system is that the 
gradient in the direction of q is small compared to the curvature associated 
with all the quantum fluctuations. In this case, an adiabatic limit results in 
which relaxation of the unconstrained meron pair slowly evolves through a se- 
quence of quasi-stationary solutions, each of which is close to a corresponding 
stationary constrained solution. Detailed comparison of the quasi-stationary 
and constrained solutions shows this adiabatic limit is well satisfied. There- 
fore, the action of a meron pair as it freely coalesces into an instanton is 
presented here. To compare with the patched Ansatz, the meron separation 
d and radii r and R of the lattice configuration are required. These are found 
by first measuring the separation of the two maxima in the action density 
A = |Mi — Miiil (using cubic splines) and their values (which are the same in 
the symmetric case, denoted by S'max = ^S/g^w'^), and then using Eq. (||) to 
give 

In Fig. ^, the total action of a meron pair is plotted as a function of d/r 
(solid line), which for the analytic case is given by Eq. (Q). Also shown are 
cooling trajectories for four lattice configurations with different initial meron 
pair separations. Each case starts with a patched Ansatz of a given separation 
and the lattice action agrees with that of the analytic Ansatz. 

The first few cooling steps primarily decrease the action without changing 
the collective variable (which is now effectively d/r), as observed above in 
Fig. ^. Further cooling gradually decreases the collective variable, tracing out 
a new logarithmic curve for the total action (dashed line in Fig. ^), which 
is about 0.25 smaller than the analytic case. This curve represents the total 
action of the smooth adiabatic or constrained lattice meron pair solutions. 
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Figure 4. Action of a regulated meron pair as a function of d/r for initial lattice configu- 
rations of oiinit = 12, 14, 18, 22 and ri^jt = 4 on a 32'' X 40 lattice. The circles show cooling 
trajectories with each cooling step marked. 



The essential property of merons that could allow them to dominate the 
path integral and confine color charge is the logarithmic interaction Eqs. (^,0) 
which is weak enough to be dominated by the meron entropy. EJ Hence, the key 
physical result from Fig. ^ is the fact that smooth adiabatic or constrained 
lattice meron pair solutions clearly exhibit this logarithmic behavior. 

In summary, stationary meron pair solutions have been found on the 
lattice exhibiting a logarithmic interaction and fermion zero mode localized 
about the individual merons in agreement with analytic results. The presence 
and role of merons in numerical evaluations of the QCD path integral should 
be investigated, and the study of fermionic zero modes and cooling of gauge 
configurations are possible tools to do so. 
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